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Context. Many haloes of nearby disc galaxies contain faint and extended features, including loops, which are often interpreted as 
i-Q ■ relics of satellite infall in the main galaxy's potential well. In most cases, however, the residual nucleus of the satellite is not seen, 

^ 1 although it is predicted by numerical simulations. 

Aims. We test whether such faint and extended features can be associated to gas-rich, major mergers, which may also lead to disc 
rebuilding and thus be a comer stone for the formation of spiral galaxies. Our goal is to test whether the major merger scenario can 

■ provide a good model for a particularly difficult case, that of NGC 5907, and to compare to the scenario of a satellite infall. 
Methods. Using the TreeSPH code GADGET-2, we model the formation of an almost bulge-less galaxy similar to NGC 5907 (B/T 
< 0.2) after a gas-rich major merger. First, we trace tidal tail particles captured by the galaxy gravitational potential to verify whether 
they can form loops similar to those discovered in the galactic haloes. 

Results. We indeed find that 3:1 major mergers can form features similar to the loops found in many galactic haloes, including in 
NGC 5907, and can reproduce an extended thin disc, a bulge, as well as the pronounced warp of the gaseous disc. Relatively small 
bulge fractions can be reproduced by a large gas fraction in the progenitors, as well as appropriate orbital parameters. 
Conclusions. Even though it remains difficult to fully cover the large volume of free parameters, the present modelling of the loops 
in NGC 5907 proves that they could well be the result of a major merger. It has many advantages over the satellite infall scenario; 
e.g., it solves the problem of the visibility of the satellite remnant, and it may explain some additional features in the NGC 5907 halo, 
as well as some gas properties of this system. For orbital parameters derived from cosmological simulations, the loops in NGC 5907 
can be reproduced by major mergers (3:1 to 5:1) and possibly by intermediate mergers (5:1 to 12:1). The major merger scenario thus 
challenges the minor merger one and could explain many properties that haloes of spiral galaxies have in common, including their red 
colours and the presence of faint extended features. 

Key words. Galaxies: evolution - Galaxies: spiral - Galaxies: individual: (NGC 5907) - Galaxies: interactions 
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!^ . 1 . Introduction vive a major merge r and form a new disc (IB antes & Hernquisl 

t4- ■ EMI 

arnesl l2002). By observing and an alysing the progenitors 
• ' The hierarchical galaxy formation theory assumes that galax- of local spi rals at intermediate redshift, Ham mer et all d2005l 
J . ies have been assembled through mergers of smaller entities. I2007I I2009h suggest that many present-day discs could have 
^ " The relative impact of b oth minor and major mergers is s t ill be- been rebuilt after a major merger. Recent theoretical advances 
I \ ing intensively debat ed (Fakhouri & Ma, 2008: [Lin et al.L I20041 also support this scenario. Many hydrodynamical simulations 
• •■ 2008 I ; [Hopkins et ail I20 08, 2009t iBell et aU l2006t lLotz et al.L indicate that new discs can be rebuilt after a major merger 
m ^ \ 120081 iBridge et all 120071) Peculiar galaxies are mostly related between disc galaxies (iRobertson et all l2006t iHopkins et all 
k>( . to merger or interaction events. For example most luminous in- 120091; [Springel & Hernquisttl2005l) . provided the initial gas frac- 

■ frared galaxies (LIRGs; Sanders & Mirabel, 1996) are merging tion is sufficiently high. In this process gas plays the key role 
d \ or interacting galaxies, as re vealed by their detailed substruc- in the disc rebuilding process. High redshift galaxies have much 

tures with multiple nuclei (see lHaan et al.ll2TJTTl) . lToomre1(ll977l) higher gas fraction than local spirals, consistently with the above 
first pointed out that a major merger could result in the forma- scenario. 

tion of elliptical galaxies. Later simulations showed that major T . ■ . , . , . ■ , f e , , , • ■ 

/ & . .,, , , . r It is expected that about half of the local disc galaxies have 

mergers of gas-poor progenitors will destroy the disc and form a , ■ ■ . n „ jzr — r — — n 

alax with a lar e s heroidal com onent exper ienced a major merger in the past 9 Gyr (IPuech et al.l 

ga axy wit a arge sp eroi a component. |201 ll) . and it should thus be expected to find imprints of such 

Partly due to their fr equency, even minor mergers chal- events within the remnant halo. For example M31, the largest 

lenge the survival of discs (Toth & Ostrikerlll992l : IStewart et all spiral of the Local Grou p, contains many faint structures in its 

2009). Thus the number of unperturbed discs at low redshift ar- halo. The Giant Stream (Ibata et al., 2001) is the most famous 

gues that there should be some mechanisms that either preserve one, together with the large north-east loo p with metal-poor 

the disc from being desttoyed or allow it to be rebuilt after the stars r ecently discovered by the PANDA team ( Richard son et all 

merger event. Analyses of merger remnants show that discs can 1201 ll) . Most of these faint features in the halo of M31 have often 

survive from a merger. Simulations with low gas fractions (< been thought to be caused by numerous minor mer ger events, 

10%) also show that a significant fraction of the gas can sur- as many as the stream number. [Hammer et al . (2010), however, 
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show that most of these faint features could, instead, be related 
to a single major merger event. In their models, most of the faint 
structures, including the Giant Stream, are due to stars coming 
back from tidal tails, a process that can be maintained for several 
Gyr. 

If the disc of many spirals is indeed rebuilt after a major 
merger, it is expected that tidal tails can be a fossil record and 
that ther e should be many loops and streams in their haloes. 
Recently Martmez-Delgado et al.1 (1201 Ol) have conducted a pilot 
survey of isolated spiral galaxies in the Local Volume up to a low 
surface brightness sensitivity of ~ 28.5 mag/arcsec 2 in V band. 
They find that many of these galaxies have loops or streams of 
various shapes and interpret these structures as evidence of mi- 
nor merger or satellite infall. However, if these loops are caused 
by minor mergers, the residual of the satellite core should be de- 
tected according to numerical simulations. Why is it hardly ever 
detected ? 

The above question is the starting point for this paper, which 
intends to test whether a major merger scenario could, or could 
not, reproduce the observations of faint loops in nearby galaxy 
haloes better than the minor merger one. We choose to study 
here the NGC5907 galaxy and its faint loop system observed by 
iMartfnez-Delgado et al.l d2008l hereafter M08) because it is of- 
ten regarded as the best evidence of a minor satellite interaction 
with a late-type spiral, i.e. not an easy configuration to be repro- 
duced by a major merger. 

The paper is organised as follows. In Sect. [2] we describe 
the properties of NGC 5907 and Sect.[3]summarises the peculiar 
features found in its halo. Simulation methods and initial condi- 
tions are described in Section|4] and in Section[5]we present the 
results of our simulations of both the galaxy and the associated 
loops. In Section [6] we discuss whether a major merger can re- 
produce a galaxy with similar properties to NGC 5907, and we 
then compare the relative merits of major and minor mergers. In 
appendices, we present the implementation of star formation in 
the GADGET-2 code (Appendix A) and we give a general de- 
scription of the loop properties formed during a major merger 
(Appendix B). 



2. Properties of the NGC 5907 galaxy 

NGC 5907 is a nearby Sc type spiral. There is a large uncer- 
tainty in the distance of NGC 5907 due to its peculiar mo- 
tion. Zepf et al. (2000) got a distance 13.5 +2.1 Mpc from 
a combination of the Tully-Fisher relation in both optical and 
near infrared and a peculiar motion model, and they adopted a 
round value 14 Mpc, which is used by M08. For consistency, 
we use the same distance (14 Mpc) as M08. The measured disc 
scale-length of NGC 5907 is 61 kpc i n #-band, while it is 3.82 
kpc in //-band dMiller & Rubinl ll995l) when scaled to 14 Mpc. 
ISaha et al l d2009l) got 3.86 kpc by fitting S pitzer IRAC 4.5 fim 
data. Bar nabv & Thronson I ([1992) obtained in H band a scale- 
height about 0.41 kpc and a bulge-to-disc luminosity ratio (B/D) 
of 0.05. 

There are still considerable uncertainties in the stellar mass 
determination of NGC 5907. We need to re-estimate it care- 
fully, accounting properly for the effec ts linked to the IMF (see 
iPuechet all 12008: Hammer et al. , 2009) Our stellar mass deter- 
mination follows the method of iBell et al ] (120031) using the K- 
band observation fr om 2MASS and an optical colour, (B - V) 
from lJust et al.l (120061) . The K s apparent magnitude is 6.757, and 
thefl-V colour is 0.86. 



So the total stellar mass is 



log l0 M+/M o = -O.4x(M^-M /fs0 )+(-O.2O6)+O.135x(B-y)(l) 



where Mk s q is the solar absolute magnitude dBell et all l200l . 
Then the total stellar mass of NGC 5907 is 6.57 x 10 10 M G as- 
suming a "diet" Salpeter IMF. The gas mass of NGC 5907 i s 
about 1.94 x 10 10 M including HI, H 2 and He (lJust et alll2006T) . 
So the total baryonic mass is 8.51 xlO 10 M , and the gas fraction 
is 23%. 

The large scatter of values found in the literature are 
mainly due to the different methods used to calculate the stel- 
lar mass. For example, methods assuming th at NGC 5907 lies 
just on the baryon ic Tully Fisher relation (e. g. iMcGaugh L 12005 ; 
Stark et all 120091) are not sufficiently accurate for an individual 
object. On the other hand, our estimate of the stellar mass may 
underestimate the contribution of the bulge because NGC 5907 
is seen edge-on and even at near-IR wavelengths, dust lying in 
the disc may considerably affect its luminosity, and its precise 
contribution to the mass is still uncertain. 



3. Peculiar properties of the NGC 5907 halo 

Even though NGC 5907 is a member of the 396th Lyon Group 
of Galaxies (LGG 396), all of the identified group members are 
at ver y large separation and co uld not be interacting with NGC 
5907 dlrwin & Maddenl 120061) . This galaxy may be considered 
as a typical isolated and an almost bulge-less spiral in the local 
universe. 

However, recent observations show that the NGC 5907 halo 
has many peculiar features. 

- There are t wo giant loops lying near the polar disc plane of 
NGC 5907. IShang et al.l dl998h first discovered one half of 
one loop, and more deep observations by M08 confirmed 
this result, revealing two giant loops. The two loops ex- 
tend up to 50 kpc from the galaxy centre. This strange 
structure indicates that NGC 5907 has experienced a former 
merger/interacting event, which could be relatively recent 
according to M08. The surface brightness of these loops is 
about 26.8 mag/arc-second 2 , corresponding to a stellar mass 
of approximately 3.5 xlO 8 M with a stellar mass surface 
density of 0.32 M pc~ 2 . 

- The gaseous and the stellar discs are strongly warped. 
Previous observations have shown that the disc of NGC 
5907 is sharply truncated at ~ 24 kpc. The radius profile 
clea rly shows a break, w hich is common to many local spi- 
rals (Ivan der Kruitt 12007). The deep observations of M08 re- 
veal that the star light extends out to the nomina l cut-off ra- 
dius. The stellar disc is also warped (see also IShang et all 
1998) in the same direction as the HI gas warp, suggesting a 
common origin for both the stellar and gas w arps. 

- The h alo of NGC 5907 is red according to iLequeux et al.l 
(1998) with V-I=1. 4 and B-V= l at 5 6 kpc off the disc 
plane. According to Zib etti etlril ([2004, see their Fig. 9) this 
is slightly redder than the stacked value of about a thou- 
sand SDSS galaxies. This implies that NGC 5907 shares the 
red halo proper ties of many nearby spiral galaxies, includin g 
M31 (see e.g. iMouhcine et all 120051: iHammer et all 120071) . 
indicating that its halo includes a significant fraction of metal 
enriched stars. 

- Mid-infrared observations revealed the presence of consider- 
able material on both sides of the disc, up to 10 kpc from the 
galaxy centre, including PAH emission dlrwin & Maddenl 
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2006) mole cular gas (| L aine et all 1201 Oh and diffuse dust 
emission (lBurgdorfetalll2009l) . This could indicate that the 
NGC 5907 disc is not completely relaxed. 

Taken together, these abnormal features suggest that a former 
merger occurred in NGC 5907. Several models were set up to 
explain some of these features. 

There are a few existing models for NGC 5907, which have 
set up to explain the formation process of the red halo or of 
the loops. To keep the disc undisturbed, a min or merger with 
a high mass ratio has generally been proposed. iLequeux et al.l 
(1998) model the red halo formation, assuming that it is from a 
merger with a red dwarf ellip tical with mass of a few 10 9 M . 
Reshetnikov & Sotnikova (2000) model the half ring discovered 
bv lShang et al.l d!998l) by using a model with a 1000:1 mass ra- 
tio. A more recent model by M08 can reproduce both loop struc- 
tures, assuming that they formed through a single merger event 
with a mass ratio of 4000:1. The small satellite, a spheroid, is 
assumed to fall into the potential well of NGC 5907, forming a 
tidal tail that literally draws out the two loops. Even though their 
model reproduces the geometry of the two loops well, it is un- 
clear whether it can reproduce all the abnormal features of NGC 
5907 described above. It is also confronted with the problem of 
the visibility of the progenitor residual. 

We therefore investigate whether a gas-rich, major merger 
can reproduce most of the abnormal features observed in NGC 
5907. Such an event may dramatically affect the whole struc- 
ture of the main galaxy, and its modelling should also reproduce 
the massive, inner components (bulge, disc) of NGC 5907. To 
achieve this goal, we need to consider a much larger number 
of constraints than that for a very minor event. For comparison, 
we list in Table Q] the constraints and their possibly associated 
parameters in both major and minor merger cases. It naturally 
leads to a huge parameter space to be investigated for modelling 
NGC 5907 by a major merger, and a natural difficulty of get- 
ting an accurate model. Our goal is to show whether a single 
gas-rich merger event can reproduce the global structure of the 
NGC 5907 bulge-less galaxy, as well as its exceptional system 
of loops in its halo. 



4. Simulations and initial conditions 

We use the publ icly available version of the GADGET-2 code 
(ISpringeli 120051) . which is a parallel TreeSP H code employ- 
ing th e fully conservative formulation (Spring el & Hern auist. 
2002) of smoothed particle hydrodynamics (SPH). In this code 
both energy and entropy are cons erved even when the s moofh - 
ing lengths evolve adaptively (Springel & Hernauist, 2003). 
Implementation of star forma t ion, c ooling, and feedback has 
been done follo wing Cox et al. ( 20061) . and our code is very sim- 
ilar to the one of lCoxetalJ (12006) (see Appendix A). 

We begin our simulations with a total baryon mass equal to 9 
x 10 10 M . This value is a few percent higher than the total bary- 
onic mass of NGC 5907, to ac count for the mass loss during the 
" 2010h. The initial conditions are set up 
(Eoioh . Each progenitor is composed 



merger ( Hammer et al 



following Ham mer et al.1 
of a stellar and a gas disc embedded in a dark matter halo. Both 
discs have an expo nential distribu tion. The dark matter model 
is chosen following [Barnes! (12002). The density profiles of these 
two components are 



Pdisc °c exp(-R/R disc ) sech (z/z d isc) . 
Phaio °c (r + a h ai o r 4 - 



(2) 
(3) 



Following lCox et al . (2006), the gas-disc scalelength is set to 
three times that of stellar disc, since observations show that the 
gas disc has a more extended distribution than the stellar disc 
(van der Kruitl 120071) . The initial stellar disc scalelength is cho- 
sen following the sc aling relation betwee n rotation curve peak 
and disc scalelength (Hamm er et al.L[2007b . We notice that fitting 
the rotation curve leaves a considerable degeneracy in the dark 
matter mass fraction in a galaxy such as NGC 5907. Here we 
adopt nine percent of the bary onic matter, which is slightly more 
than other adopted values dHoekstra et all 20051: iDutton et all 
2010) but less than those a dopte d in Barnesl d200ST ICox et al.l 
(2006]), and lHammer et al l (|201Q|) . The final scale length of the 
dark matter profile is fix ed using the baryonic Tully-Fisher rela- 
tion dPuech et al.Ll2.010b . The disc scaleheight is s et to be ten per- 
cent of the scalelength. We follow the method o f B arne i d200l 
to build the initial galaxy, which is stable enough for current 
work (see Appendix A). Parameters for each model are listed in 
Tabled 

The adopted mass distrib ution of halo is more susceptible 
to tid al tail formation dDubinski et all Il999t iD'Onghia et all 
2009). To test the effect of the h alo profile on the final result, 
we considered a Hernquist model (Hernauisl. 119901). This model 
is clo se to an NFW profile in the inner region ( Springel et all 
2005) and we find that the results do not change significantly 
(M3L34H model in Table 0. The halo concentrations for the 
primary and secondary interlopers are found to be C=l 1.48 and 
15.6, respectivel y, the first value being very close to that of 
ICox et all d2006l) . who model a gas-rich Sbc galaxy. We did not 
assum e halo rotation since it is not expected (Springel & White, 
1999) to have a strong influence on tail formation. Simulations 
use between 180000 and six million particles to test the effect of 
resolution. 

The gas fraction is assumed to be quite high (from 60 to 
80%), for progenitors that are intermediate mass g alaxies at z> 
1. For comparison gas fractions of 60% are found bv lDaddi et al.l 
(2010) for z~ 1.5 galaxies with stellar masses of 5xlO lo M and 
similar values are found for massive galaxies at z~ 2 (Erb et al, 
2006; see also Rodrigues et al. 2012 in preparation). To evaluate 
the mass ratio betwee n the two merging galaxies, we follow the 
Hopk ins et al.l d2010h zeroth order scaling relation between gas 
fraction, mass ratio, and bulge fraction. The small bulge compo- 
nent of NGC 5907 gives a constraint on the mass ratio, which 
cannot be too small and should be larger than or equal to 2. 

By construction, the present simulations assume a single en- 
counter of two galaxies, without accounting for other external 
supplies of gas such as cold flows. Not being able to extract a 
galaxy similar to NGC5907 from cosmological simulations has 
several consequences. For example, the assumed gas fraction in 
the progenitors is certainly an upper limit, because some addi- 
tional gas may be captured from the IGM before fusion. The sta- 
bility of the gas-rich progenitors when they are isolated has been 
verified in Appendix A, and do not account for gas accretion 
and clumpy fragmentation as could be inevitable with sufficient 
levels of gas accretion. Besides this, our simulations of an iso- 
lated disc reproduce those of Cox et al. for all configuration of 
density and feedback. Possibly the absence of clump is related 
to ours and the Cox et al. recipes for SPH simulations ; how- 
ever, recent results from the AREPO code dKeres et all 1201 lh 
show galaxies without numerous dense gaseous clumps, possi- 
bly contradicting with the numerous clu mpy, intermediate-mass 
galaxies at z~ 2 dElmegreen et all 120091) . On the other hand, it 
is unclear whether clumpy disc progenitors would have changed 
our results after a major merger, which mostly destroy progeni- 
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Table 1. Comparing the number of constraints and parameters between minor and major mergers. 



Constraints 


Major merger 


Minor merger 


Model parameters 


Loops shape 


Y 


Y 


inclination, orbital parameters, view angle 


Loops mass 


Y 


Y 


mass ratio, gas fraction 


B/T 


Y 


N 


mass ratio,star formation history, gas fraction 


Rotation Curve 


Y 


N 


initial conditions 


Gas warp 


Y 


N 


inclination, orbital parameters 


Disk scalelength 


Y 


N 


inclination,orbit parameters, mass ratio 



tor's features, except possible bulge-like structures th at could be 
enhanced in clumpy galaxies dBournaud et al.L [201 lb . 

For some of our models, we assume a star formation history 
with a varying global efficiency in transforming gas to stars, in 
order to preserve enough gas from being consumed before fu- 
sion. Although this fine-tuned star fo rmation history may have 
some physical motivations (see e.g. Ham mer et all 1201 Ol) . its 
main role is also to ensure the formation of stars after the emer- 
gence of the gaseous disc just after fusion. 

We adopt a retrograde-prograde mer ger, which is favourabl e 
for building a disc after the merger dHopkins et al 1 l2009h . 
The rotation direction for the secondary progenitor is set to 
be prograde and that of the massive progenitor, retrograde. 
This choice is because resonances have a strong effect on 
long tidal tail formation, which is e nhanced in prograde and 
suppressed in retrograde e ncounters |Toomre & Toomrel 1 1972b 
Barn es & Hernquistt 119921) . We also find that in order for the 
loop shapes to match the observations, the orbital eccentricities 
should not be as low as found by M08, who used an orbital ec- 
centricity of 0.42. We find, however that an orbit eccentricity of 
0.9 is enough to have loop sha pes that match the observatio ns. 

Cosmological simulations dKhochfar & Burk ert. 2006) pro- 
vide a relationship between the ratio of the pericentre to the 
virial radius and the orbit eccentricity. We carefully choose peri- 
centre and eccentricity so that our orbital geometry fulfils this 
relation. For example, we estimate the virial radius of the pri- 
mary galaxy halo at the beginning of the simulation (8-9 Gyr 
ago, z>l) to be about 130 kpc from the relationship between 
virial r adius and total mass der ived from cosmological simu- 
lations (Springel&_WhM[l999). Assuming a pericentre of 25 
kpc implies eccentricity v alues of above 0.88 (see figure 6 of 
iKhochfar & Burkertl 120061) . 

5. Modelling NGC5907 and its loop system 

5.1. Formation of giant loops 

5.1 .1 . Mechanisms of loop formation during the merger: 
setting the orbital parameters 

Appendix B describes the basic properties of loops formed in a 
major merger from tidal tail particles, which later on are captured 
by the gravitational potential of the galaxy. Here we use these 
generic properties of tidal tails and loops to investigate different 
viewing angles and different forming epochs to reproduce the 
loops found by M08 in the NGC 5907 halo. 

The assumption that the loops are well fitted by ellipses pro- 
vides important constraints on the viewing angles. To form long- 
lived loops, their associated particles should not interact with 
the newly formed disc. We notice that one loop (noted as SW- 
SE-E1-W2 in M08) intersects the optical disc. Assuming it is a 
projection effect, this implies that the angle between the line of 
sight and the loop plane should be 60 degrees or larger. This is 
consistent with the M08 model, which assumed a loop rotation 



by 57 degrees. Moreover the inclination of the loop plane to the 
disc plane is constrained by the observed tip position of the el- 
lipse. The observed angle between the tip of the SW-SE-E1-W2 
and the direction perpendicular to the disc, is quite small (less 
than 15 degrees), indicating that the loop plane should be almost 
polar to the disc. The loop shape also depends on the precession, 
intrinsic eccentricity, and position. 

Table[2]summarises the properties of the models that are used 
in this paper, including initial parameters, resolution, star forma- 
tion, or feedback history, as well as the structural parameters of 
the remnant galaxy well after t he merger. Feedba ck model 1 uses 
five times median feedback of ICox et alJ J2006) used. Feedback 
model 2 (M3L34F) uses five times this median feedback before 
fusion,and changes to a low feedback after fusion. The bulge 
Sersic index n and disc scalelength are obtained by fitting the 
surface brightness density distribution with a Sersic profile, and 
two exponential functions. Model M3L34G6 changes gas frac- 
tion. Model M3L34H uses Hernquist model. Model M3L34A, 
M3L34B are similar to M3L34 after a fine tuning of galaxy in- 
clinations. Figure [TJ gives an example of a 3:1 merger and dis- 
plays its evolution at different epochs. The first passage occurs at 
about 1 .0 Gyr after the beginning of the simulation, and second 
at 2.7 Gyr. The first tidal tail is formed after the first passage, 
and expands at 2.7 Gyr, i.e. just before the second passage, as 
shown by a yellow arrow in the upper third panel. We found that 
by increasing the pericentre radius, this tidal tail becomes less 
prominent and provides fewer stellar particles falling back to the 
newly formed galaxy. After the second passage, a second tidal 
tail is formed, as shown by a black arrow in the upper right hand 
panel of Fig. Q] (3.5 Gyr). Particles in this second tidal tail may 
form loops, and at ~ 4.2 Gyr, the first loop begins to form (red 
arrow in the second left row panel). At this time the loop size is 
small. As time elapses an increasing number of high-elevation 
particles come back from the tail, and the first loop expands (see 
Appendix B). At ~ 5.6 Gyr, the new disc is well formed, and 
a second loop appears (pink arrow in the 2nd row 2nd column 
panel). As time evolves, more and more loops appear. At 8.6 
Gyr, the third and fourth loops (green and blue arrows) are well 
formed, and their structures show some similarities with those 
of NGC5907. At the same time, the first and the second loops 
have considerably expanded, which unavoidably leads to their 
dilution. 

The two observed loops have approximately the same appar- 
ent size, which we interpret to be a projection effect. Indeed, the 
higher order loop should be smaller than the lower order one, 
as described in Appendix B. A colour bar in the top left hand 
panel of Fig. Q] indicates the stellar-mass surface density. The 
surface mass density of the loops is about ~ 0.1 to 1 M pc" 2 . 
M08 roughly estimate the loops' surface brightness as ~ 0.32 
M Q pc~ 2 . Even though there is large uncertainty in the observed 
values, the model is able to reproduce the order of magnitude of 
the loop surface brightness. 
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Table 2. Eight models used in this study and their associated parameters. The parameters above the bar describe the progenitors, 
while those below it describe the merger remnant. 



parameters 


M3L34G6 


M3L34 


M3L12 


M3L23 


M3L34H 


M3L34F 


M3L34A 


M3L34B 


mass ratio 


3 


3 


3 


3 


3 


3 


3 


3 


halol core size 


11.8 


11.8 


11.8 


11.8 


20 


11.8 


11.8 


11.8 


halo2 core size 


6.5 


6.5 


6.5 


6.5 


11 


6.5 


6.5 


6.5 


stellar disci scalelength 


4.6 


4.6 


4.6 


4.6 


4.6 


4.6 


4.6 


4.6 


stellar disc2 scalelength 


3.5 


3.5 


3.5 


3.5 


3.5 


3.5 


3.5 


3.5 


Gall incx 


-147 


-150 


-130 


-140 


-150 


-150 


-145 


-155 


Gall incz 


-175 


-180 


30 


-70 


-180 


-180 


-170 


-180 


Gal2 incy 


-30 


-30 


-30 


-30 


-30 


-30 


-30 


-30 


Gal2 incz 
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Fig. 2. Inclination angle between disc, loop, and orbital planes 
as a function of time. The lower solid lines show the angle be- 
tween the loop and the orbital planes. Lines within shaded re- 
gions indicate the inclination angles of the disc plane relative 
to the orbital plane, while the line-shaded region around 90 de- 
grees indicates the angles between disc and loop planes that are 
consistent with the observations. Different colour lines show the 
simulations with different initial disc inclinations to the orbital 
plane. All the models used here have a mass ratio of 3 (from 
bottom to top, models M3L12, M3L23, M3L34B). 



Figure [2] shows how the disc and loop inclinations evolve 
with time. The loops are formed by particles coming back from 
the tidal tail. We quantified the loop plane inclination relative 
to the orbital plane from the angular momentum of the tidal 
tail particles. Because the angular momentum of the tidal tails 
is dominated by the orbital angular momentum, tidal tail par- 
ticles mostly lie on the orbital plane, and at most 20 degrees 
away from it. To quantify the disc plane, we use the angular mo- 
mentum of the post -fusion stellar particl es, following the pro- 
cedure presented in iHammer et all (l2010h . With evolving time, 
the disc plane systematically rotates so as to be aligned to the 



orbital plane. Using the inclinations of disc and tidal tail parti- 
cles relative to the orbital plane, we can determine the inclination 
between the loops and the disc plane that is consistent with ob- 
servations. The angle between loop and disc plan es is about 80 
degrees in M08, while iReshetnikov & Sotnikoval d2000l) use an 
angle of 90 degrees to fit one loop. We considered this ten de- 
gree difference as the uncertainty in the angle between the loop 
and disc planes, as well as the uncertainty in our models, so in 
Fig. [2] the acceptable angles are shown by the line-shaded re- 
gion. Because higher order loops appear at later epochs, the ini- 
tial disc inclination angle should be larger to make the disc and 
loop plane close to polar at later times. 

5.1.2. Comparison with observations: which loops matter ? 

As shown in the previous section, a major merger can gener- 
ate loops with a geometry similar to the observed ones, at a 
time that depends on the initial orientation of the main progeni- 
tor disc plane. Because observations are surface brightness lim- 
ited, we illustrate in Fig. [3] how the visual impression depends 
on the surface-brightness threshold. By varying the threshold, 
the thin disc, thick disc, and loop features are progressively re- 
vealed. Colour bars in Fig. [3] indicate that loop regions with sur- 
face brightness ranging from 0.1 to 1 M G pc~ 2 coincide quite 
well with the observed loops. Since the loop surface brightness 
is 0.32 M pc~ 2 , this argues that observations identify only the 
brightest part of the loop system, i.e. mostly above the galaxy in 
the panel corresponding to 8.6 Gyr. 

After the merger, particles coming back from the tidal tail 
continuously form loops of increasing order as illustrated in Fig. 
[TJ During this process, while higher order loops are forming 
and reach sizes comparable to the observed loops, lower order 
loops continuously expand and become more and more diluted 
with time. Another factor favouring the low-order loop dilution 
is related to the tidal tail surface brightness distribution, which is 
fainter farther from the galaxy. The low-order loops are formed 
at any time by particles coming back later, which results in lower 
surface brightness compared to high-order loops. 

We also find that a more inclined initial disc can capture the 
material of the secondary interloper better by making it rotate 
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Fig. 1. Star and gas particles distributions for a merger with mass ration 3 (model M3L34A) at different epochs. Star particles are 
shown in the top two rows. The fit to best observations is at 8.6 Gyr, at which the structures of third and fourth loops can match the 
observations. The gas surface density is shown in the bottom two panels with gas density indicated by the colour bar. The size of 
each panel is 300 by 300 kpc and the arrows are explained in the text. 



closer to the disc plane. Thus, there is little residual signature 
when the disc is seen in an edge-on position and favours ob- 
served loops of a high-order type. This is illustrated well in Fig. 
|H which compares the formation of different order loops at dif- 
ferent epochs. Residuals from the secondary decrease dramati- 
cally with time, favouring high-order loops that match the ob- 
servations. Figure|4]also illustrates that lower order loops almost 
vanish with time. 

In our model, the left loop has a higher order than the right 
one. This is in good agreement with observations, since the left 
loop has been discovered first bv lShang et alJ (Il998h . and the dis- 
covery of the two loops by M08 has been done with much deeper 
observations, implying that the left loop is brighter than the right 



one and that surface brightness threshold has to be accounted for 
when modelling the loops. 

Another strong argument favouring observed loops of the 
high-order kind is related to the nature of the progenitors, which 
are assumed to be gas-rich galaxies to allow the formation of 
a significant disc. Such gas-rich galaxies are common in the 
high-redshift universe, and far less common at intermediate red- 
shifts. Using this argument, together with the absence of resid- 
uals, favours a scenario for which the observed loops coincide 
with a combination of at least third and fourth order loops, with 
a time after interaction longer than ~ 8.6 Gyr. Interestingly, this 
scenario can be falsified by further observations: deeper obser- 
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Fig. 3. Stellar mass surface density shown with different scales of stellar mass surface brightness at 8.6 Gyr for the model M3L34. 
Each panel size is 300 by 300 kpc. 
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Fig. 4. Comparing loop structures formed by combining loops of different order for mass ratios of 3. The left panel shows the first 
and second loops (model M3L23), the second panel shows the third and fourth loops (model M3L34), while the first loop appears 
at the bottom of the galaxy and becomes diluted. The size of each panel is 300 by 300 kpc. The loop stellar surface mass densities 
are consistent with observations from M08. In all these images the threshold has been chosen slightly below the observed one from 
M08. 



vations of a wider field of view should detect the fainter lower 
order loops. 

5.2. The rebuilt central galaxy: disc and bulge 

To test whether the final merger remnant is consistent with the 
disc of NGC 5907, we follow the decomposition method of 
iHammer et al ] (120101) . Each final remnant is decomposed into 
three-components, namely bulge, thick disc, and thin disc. This 
method is based on angular momentum distribution taking ad- 
vantage of the full t hree dimensional in formation for each parti- 
cle (see details in Ham mer et al.L l20loh . Figure [5] gives an ex- 
ample for a 3:1 merger model. The top panels show the an- 
gular momentum distribution of young, intermediate age, an d 
old stars, which are defined following IHammer et al.l (2010). 
Most young and intermediate age stars have angular momen- 
tum along the disc polar axis (0.9<MAz/MA <1) or the thick 
disc (0.7<MAz/MA <0.9), with the rest of them mostly concen- 
trated in the bulge. The top right hand panel shows the presence 
of many old stars at a large distance from the centre, with an- 
gular momentum very different from that of the disc. They are 
mostly particles of the loops, but not do rotate with the disc. The 



bottom panels of Fig.[5]show the projected mass distribution for 
different components, confirming the presence of a prominent 
thin disc including > 70% of the baryonic mass after the merger 
event (left panel). 

We also fitted the stellar surface brightness distribution simi- 
larly to what is usually done from observations. For each merger 
remnant, we fit the galaxy surface brightness with a bulge that 
has a free Sersic index and an exponential disc (Sersic n=l). For 
each merger remnant, we define the disc plane using the three- 
dimensional information as defined after the galaxy decomposi- 
tion using the angular momentum (see Fig . fSJ* . The stellar surface 
density is then measured in logarithmically spaced annuli in the 
disc plane. 

Figure. |6]shows one example of a fit for a model with a mass 
ratio of three. During the fit, we excluded the very central part 
because it could be affected by numerical softening. We note that 
an additional exponential component is necessary to fit the out- 
skirts. This additional component presents a surface brightness 
similar to or fainter than that of the loops and corresponds to the 
inner halo of the simulated galaxy. The results of the different 
models are listed in Table [2] The disc scalelength of our models 
ranges from 3.9 to 4.5 kpc and is consistent with the observed 
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Fig. 7. Comparison of observed and simul ated (model M3L34A) 
rotatio n curves. The blue points are from lSancisi & van Albada 
( 1987) and are based on HI observa tions. Black soli d point s are 
from Ha and Nil observations by iMiller & Rubinl (1 19951) . in- 
cluding near the central regions. Solid line shows the prediction 

from the mass model V c = J GM ^ <r \ red points are direct mea- 
surements from the gas kinematics within ±1.5 kpc from the disc 
plane. 

scalelength value of 3.84 kpc. The bulge fraction ranges from 
15 to 23%, and depends on both the model and the star forma- 
tion history. We caution that this bulge fraction is derived from 
the stellar mass distribution, while in some models ~ 40% of 
baryons are still made of gas. 

Observations show that NGC 5907 is relatively gas-rich (gas 
fraction of 23%), which can be compared to the final gas frac- 
tions of the different models listed in Table [2] Models with con- 
stant feedback result in excessive gas fraction in the final rem- 
nants. We verified that the final gas fraction may reach the ob- 
served value using a tuning of the feedback (or an increase in the 
star formation efficiency) just after th e fusion (model M3L34F), 
as assumed by Ha mmer et al.l (l2010l) . This also leads to a lower 
bulge fraction, with B/T=15%. 

The surface-brightness fitting method also gives the bulge 
half light radius and the bulge sersic index. Most of our real- 
izations of gas rich major mergers result in low Sersic index 
values (< 1.5), i.e. com parable to those of pseudo bulges (e.g., 
iKormendv & Kennicuttll2004l) . 

5.3. Rotation curve 

We also compared the model rotation curve with observations. 
We collected two observations for rotation curve of this galaxy 
as shown in Fig. [7] Even though the observations show some 
discrepancies between them, our model reproduces the observed 
rotation curve quite well. 

5.4. Disk warp and disturbed material above the disc plane 

The merger remnant is not very relaxed at the time of observa- 
tion, which makes the disc dynamically hot. This violent process 
naturally results in a disc warp and disturbances. The bottom 
rows of Fig. Q] show the gas evolution after the merger and the 
presence of a warp at 8.6 Gyr. After the merger, gas is decou- 
pled from stars owing to their different dissipational properties, 



which is consistent with the observational result that the loops 
only include stars. The merger process makes the gas warped at 
the edge of t he dis c, and this can be maintained for several Gyr. 
Shang etaUd 19981) did a 6 hr radio observation with the VLA. In 
their figure 2 of an integrated 21 cm intensity map, the gas dis- 
tribution is very warped and extends above and below the plane 
of the disc. 

The gas above the disc plane shows non-relaxed motions 
that are also exp e cted w ith our modelling of a major merger. 
Ilrwin & Madden! ((2006) analysed the PAH distribution and 
found PAH emissions well above the disc, up to 8.3 kpc. CO 
observation of the central parts of the disc show a steeply rising 
rotation curve and non-circular molecular ga s motions, presum- 
ably due to a bar (Garcia-Burill o et al ., 1997). These observa- 
tions could be consistent with a major merger origin of NGC 
5907. 

5.5. Effects related to resolution 

To test the effect of resolution, we compared simulations with 
different numbers of particles, from 180 k particles to six million 
for model M3L23. This test is used to verify whether weighing 
dark matter particles and gas particles affects our final result dif- 
ferently. In this study, the mass of the individual particles for 
dark matter, gas, and stars has ratios from 10:1:2 to 4:1:1. By 
varying particles mass ratio and increasing total particles number 
up to six million, we find that our results, including the structure 
of the loops, are not significantly affected by resolution. 

6. Discussion and conclusions 

6.1. Can the properties of NGC5907 be reproduced by an 
ancient and gas-rich major merger? 

We modelled NGC 5907 as the result of a gas-rich major merger, 
using progenitors with a mass ratio of three. With gas fractions in 
the progenitors ranging f rom 60 to 80% and a s tar formation his- 
tory similar to the one in Hammer et al .(2010), our modelling is 
able to reproduce the main observed NGC 5907 structural pa- 
rameters (the thin disc and its rotation curve, and the bulge). 
The exceptional features of NGC 5907 can be reproduced, to- 
gether with the central galaxy properties, especially if we com- 
pare the observed loops to the high-order loops expected in a 
major merger model. Given the extremely large number of pa- 
rameters, as well as the very numerous constraints provided by 
the observations, we cannot claim that we have already identified 
the exact and unique model of NGC 5907 and its halo properties. 
We nevertheless succeeded in reproducing the loop geometry, 
and a disc-dominated, almost bulge-less galaxy. The limitations 
of our modelling are two-fold, with one related to observational 
uncertainties, the other to our modelling procedure. 

Distance uncertainties of ~ 30% translate into uncertainties 
in the total baryonic mass of the galaxy, which may affect the 
accuracy of our modelling. However, the main uncertainty when 
comparing observations to modelling is related to the striking 
difference of the NGC 5907 disc at visible and near-IR wave- 
lengths as illustrated in Fig. [8] At visible wavelengths the thin 
disc axis ratio is almost eight, the bulge shows a negligible con- 
tribution, while this exposure is deep enough to detect the very 
low surface brightness loops. This is in sharp contrast to obser- 
vations at near-IR wavelengths, which reveal a thicker disc, with 
an axis ratio of 4.2, and do not detect the loops. 

Two effects can explain the above behaviour of the edge-on 
disc of NGC 5907. Firstly, it has been shown that thick discs 
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Fig. 5. The top panel gives the distribution of z component of the angular momentum of stellar particles for model M3L34. The 
galaxy has been projected edge-on using the angular momentum of the baryons. From left to right, star of different age (from young, 
intermediate, and old) are shown. Old stars (red dots) already exist in the progenitors at the beginning of the simulation. Young 
stars (blue dots) are defined with ages over 3 Gyr after the beginning of the simulation. Ages between young and old are defined as 
intermediate stars (green dots), which are formed before the fusion. The bottom panel shows the projected mass distribution of the 
thin disc (left), and the bulge and thick disc (right). The thin disc is selected from its angular momentum ratio (MAz/MA> 0.9, see 
Hammer et al, 2010), which cuts the warp off. Each panel has a size of 300 by 300 kpc. 



Table 3. Comparing major and minor merger hypotheses, where "Y" means consistency with the observations. 



Observational features 


Major 


Intermediate 


Minor Comments 




merger 


merger 


merger 




3:1-5:1 


5:1-12:1 


> 12:1 



Loop shape 


Y? 


Y? 


Y 


In minor merger, the loops trace the progenitor orbit, while in a 

major merger, loops are formed by particles coming back from a tidal tail 


Loop size 


Y 


Y 


Y 




Loop surface mass density 


Y 


Y 


Y 


Most models can reproduce the stellar mass surface density 


Loop eccentricity 


Y 


N? 


N 


With an orbital eccentricity of 0.9-1, a merger with mass ratio > 12: 1 
would need more than a Hubble time to reach fusion 


Visibility of 


Y 


Y 


N 


Major mergers provide a thick disc component after complete fusion of the nuclei, 


remnants 








and low dynamic friction means the remnant nucleus should be seen in a minor merge 


Gas warp 


Y 


Y 


N 


Major merger predicts a gas warp unlike a minor merger with high mass ratio 


Gas & PAH emission 


Y 


Y? 


N 


Major merger predicts residual gas and PAH emission above the disc, and 


above the disc 








this should be explained by another mechanism in a minor merger 


Colors of loops 


Y 


Y 


N? 


For a minor merger the satellite should be a relatively massive, red dwarf elliptical 



are pr ominent at near-IR wavelengths (see e.g. IComeron et all 
1201 lb . and secondly, the presence of diffuse interstellar dust 
up to 10 kpc above th e disc has been identified in NGC 5907 
(Bure dorf et all 120091) . At 5 kpc from the disc, we estimate the 
surface brightness of the disc to be 22.5 mag arcsec' 2 at 3.6/mi, 
which translates into 25.3 mag arcsec~ 2 in the AB system. At 
this location, there is a very faint emission in the image by M08 
with ~ 28.0 mag arcsec~ 2 in the R band, i.e. 28.3 mag arcsec~ 2 
in the AB system. Such a huge gap of three magnitudes between 
near-IR and visible light is probably caused by a combination of 



both effects, very red stars in the extended thick disc, and dust at- 
tenuation in the visible. It translates to / V (V) ~ 0.06x/ v (3.6yum). 
To estimate the impact of dust extinction at 5 kpc on both sides 
of the disc, we assume that the thick disc stars cannot be in- 
trinsically r edder than M stars. Fo r a non-extincted M star with 
r o =3500K. lBeniamin et ail d2007l see their Fig. 7) find / V (V) ~ 
0.15x/ v (3.6//m). Thus even if the thick disc was only populated 
by M stars, it should have been detected in the very deep ex- 
posure of M08, except if the visible light has been extincted by 
a factor 0.15/0.06=2.5, which corresponds to Ay=0.9 (Galactic 
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Fig. 6. Stellar surface mass density distribution as a function of radius for model M3L34. Open circles show the total stellar density 
together with a fit by different components (solid line). Sersic bulge and exponential disc components (indicated by different colours) 
are used to fit the stellar mass surface brightness. The fitting process excludes an area with a radius equal to 2 times the softening 
length (labelled by arrows, see the text). 




Fig. 8. Top: Deep observations of NGC 5907 in the visible 
light showing the thin disc with the gigantic loops (from M08). 
Bottom: Observations of NGC5907 with IRAC at 3.6/rni reveal- 
ing the whole extent of the disc that is twice thicker than in the 
visible light. 



extinction law for R=3.1). For illustration, by replacing M by 
K stars (with r c =4500K) would lead to an extinction factor of 
12 instead of 2.5, emphasising the importance of dust extinction 
on both side of the edge-on disc. Ideally, only a very deep ex- 
posure at near-IR wavelength can provide an image that can be 
compared to the stellar mass distribution of a model. 

Another substantial uncertainty from observations is related 
to the B/T estimate. Our model predicts B/T values around 20%, 
and a specific implementation of varying star formation history 
(or varying feedback history) reduces this value to 15%. The 
latter value is approximately twice the value we derive in fit- 
ting the luminosity profile of the 3.6 /mi image (B/T ~ 0.08), 
which should be less affected by dust than at visible wave- 
lengths. In fact our simulated B/T for the stellar mass may not 
coincide exactly with estimates from the luminosity. It has been 
argued, for example, that the IMF in elliptical galaxies is more 
dominated by low-mass stars than that of the Milky Way disc 
(van Dokkum & ConrovL l2010). Given the similarities between 
ellipticals and bulges, one may speculate that the NGC 5907 



bulge includes a larger fraction of low mass stars than the disc, 
leading to a B/T similar to the simulated one. 

Limitations of our modelling are linked to the huge parame- 
ter space to be investigated. It leads to loops that behave like the 
observed ones, with some differences in the size and geometry. 
In Fig. [9] the modelled loops also show some differences with 
the observed ones when comparing their locations with respect 
to the disc. This indicates that, while we correctly identify the 
mechanism for loop formation, we do not yet recover the best so- 
lution for NGC 5907, which would require an enormous amount 
of time for fine tuning the models. Other limitations come from 
(1) the use of GADGET-2, which, for given resolution and star 
formation, feedback, and c ooling recipes, pr ovides less well de- 
fined discs (see Fig. 11 of iKeres et all 120111) with lower angu- 
lar momentum than those predicted by the AREPO code. (2) 
The fact that our simulations are done in isolation, while other 
sources of gas supplies may be brought to the galaxy as done in 
cosmological simulations. Accounting for the above could have 
led to more realistic discs, and as it reduces the need for too high 
a gas fraction in the progenitors. 

6.2. Comparing the relative merits of minor and major 
merger models 

As mentioned above, the main drawback of the major merger 
scenario is caused by the enormous number of constraints to 
account for, leading to simulated configurations that still show 
some discrepancy with observations. This is in sharp contrast to 
a minor merger case for which the satellite orbit is only con- 
strained by tracing the loops, without any influence on the main 
galaxy. This point, however, is mostly technical and should be 
discarded when discussing the relative relevance of both scenar- 
ios. 

Table [3] summarises the exceptional features revealed by 
deep observations of NGC 5907 at visible, mid-IR, and radio 
wavelengths, and displays the merits of major, intermediate, and 
minor mergers, respectively, to reproduce them. To sample a 
wider range of mass ratios, we ran simulations of 10:1 and 15:1 
mergers assuming a pericentre radius of 25 kpc and orbits con- 
sistent with the formation of loops similar to the observed ones. 
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Fig. 9. Left: Deep observations of NGC 5907 in the visible 
light showing the thin disc with the gigantic loops (from M08). 
Middle: One simulation of the loops at T=8.0 Gyrs with the 
M3L34G6 model, in which we assume a dust extinction screen 
affecting the stellar light density by a factor of 2.5 on both sides 
of the edge-on disc (see discussion in Sect. 6.1). Right: Same 
model as middle panel shows all components with accounting 
for an extinction by a factor 12. The real extinction should be 
between the middle and right panel. 



While the 10: 1 merger occurs at T= 9 Gyr after the beginning of 
the simulation, the 15:1 merger would take more than a Hubble 
time to effectively merge. Thus we estimate that 12:1 is a limit 
for reproducing NGC 5907 with the orbital parameters expected 
from cosmological simulations. 

A possible advantage of the major merger is linked to its 
being likely to predict all these features together, while the mi- 
nor merger case can only reproduce the loops. Features like the 
warp or the gas and PAH emission above the disc can be ex- 
plained well by other mechanisms, such as additional interac- 
tions or ejection of dust particles from the thin disc by galac- 
tic winds. The presence of considerable amounts of non-relaxed 
material (dust, gas including in molecular phase) in the disc out- 
skirts, as well as the strong warp favours a scenario for which the 
NGC 5907 is not a relaxed disc, but this cannot lead to a decisive 
conclusion. 

Possibly the absence of a remnant signature, together with 
loop colour are more compelling for disting uing between the tw o 
scenarios. Using the IRAC 3.6yum image (Ash bv et all l200l . 
we were able to model the light distribution with GALFIT 
dPengetallfeolOl) using a bulge, plus disc decomposition. After 
subtracting the model from the image, we examined the presence 
of a possible residual of the progenitor of a minor merger at the 
locations where the loops intersect the edge of the disc. At such 
a large distance (20 kpc) from the disc centre, it is unlikely that 
dust can affect the light that much, so we find no object that could 
be responsible of such an event. This is really problematic for a 
minor merger scenario, because with such an extended orbit, it 
is unlikely that tidal forces could destroy it. In the major merger 
case we have shown that the low surface-brightness features left 
by the secondary interloper gradually vanish with evolving time, 
because they are progressively captured by the disc (see Fig.|4|i. 

The colour of the loops has been e stimated by differ- 
ent groups and con firmed to be both red dZheng et all [l999; 
iLequeux et all 1 19981) a nd comparable to that of elliptical galax- 
ies dZibettiet all 120041) . With the red colour R - I ~ 0.5 ± 0.3, 
IZheng et alJ (119991) found it is consisten t with Galactic globu- 
lar clusters with [Fe/H] ~ -1. Following Hammer et al. (2007), 
we show in Fig. \W\ the relation between rotational velocity 
and colour-inferred iron abundances of the stellar disc outskirts. 
NGC 5907 follows the same relation as defined by most spi- 
rals, which could be easily explained if stars in the inner halo 
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Fig. 10. Comparing iron abund ances of outskirts of N GC 5907 
with that of other galaxies from Ham mer et alJ d2007[). The iron 
abund ances of the NGC 5907 inner halo are from IZheng et alJ 
(Il999h . who estimate them by comparing the loops colour with 
that of Galactic globular clusters. 



are e nriched by an ancient major merger (see lHammer et all 
2007). Given their location in the mass metallicity relation, a 
dwarf spheroid origin for the loop progenitor would generate 
a bluer colour than what is observed. For example, the colour 
of Sagittarius and its stream [Fe/H] =-1.2 dSesaii 1201 ll) hardly 
matches the colour of the NGC5907 halo. 

Assuming a cosmological origin for the encounter that 
is responsible for the loops has important consequences. In 
fact, this implies orbit ecc entricities in excess of e=0.85 (e.g. 
iKhochfar & Burkerl 120061) for a pericentre larger than 25 kpc. 
This would exclude the M08 model that assumes e=0.42, as well 
as all models with mass ratio over 12:1, simply because the two 
galaxies would not have time to merge within a Hubble time, 
leading the problem of the absence of progenitor residual too 
critical. 

Examination of Table [3] confirms the above: only mergers 
with mass ratios between 3:1 and 12:1 are consistent with most 
of the observations of the NGC5907 halo. Given the observa- 
tional constraints and our modelling, it is not possible to distin- 
guish between the merits of intermediate mergers (5:1 to 12:1) 
and major mergers (3:1 to 5:1). In fact, we have failed to re- 
produce the observed loops with a 10:1 merger, although we 
certainly did not investigate the whole range of parameters to 
exclude this. On the other hand, with orbital eccentricities from 
0.9 to 1, less dynamical friction is expected for a higher mass 
ratio, and it could be difficult to match the low loop eccentricity. 

Possibly deeper exposures on larger scales of the NGC5907 
halo could reveal many very faint new structures that may help 
disentangle between these two alternatives. Using state-of-the- 
art AREPO code to reproduce the NGC5907 loops would also be 
useful for providing a more realistic the disc and the distribution 
of matter in its surrounding. 

6.3. Concluding remarks 

This study shows that a major merger scenario for the forma- 
tion of giant loops in NGC 5907 may explain the observations, 
as well as the infall of a minor merger with a mass ratio lower 
than 12:1. Together with the rebuilding of a thin disc in an al- 
most bulge-less galaxy, a gas-rich major merger forms gigantic 
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loops, within which stars may orbit in loops for several billion 
years. Low-mass ratios down to three are particularly efficient in 
removing the residuals of the secondary interloper, and are still 
consistent with the observational uncertainties on the bulge-to- 
total mass ratio. This is especially true when high order loops 
are considered, several billion years after the merger. This con- 
trasts with minor mergers that are assumed to orbit at large dis- 
tance from the galaxy centre, and should have left imprints of 
a residual core. This proposition of a major merger origin could 
apply as well to the faint tidal features found in many other spiral 
galaxies. Interestingly, our model of NGC 5907 can be falsified 
by new and extremely deep observations of a wider field sur- 
rounding this galaxy: if NGC 5907 has been formed by a major 
merger, there should be faint and extremely large structures at 
hundred kilo-parsecs or more, coinciding with lower order loops. 

Future work will include modelling other nearby spiral 
galaxies with large and faint, extended features in their haloes. 
Most of them show no trace of a residual core of a minor merger 
residual. While it is accepted more and more t hat bulge-less 
galax i es could resu lt from g as-rich maj o r mer gers (G uedes et al.L 
1201 U iBrook et all 1201 U iFont et all 120111) , the consecutive 
modelling of faint features in their halo should put considerable 
constraints on the orbital parameters. A wide diversity of orbital 
parameters has been already found by [Hammer et al.l (120091) in 
modelling distant, unrelaxed galaxies by major mergers. These 
distant galaxies are likely similar to the progenitors, six bil- 
lion years ago, of present-day spirals, and linking them together 
could provide another crucial test for the spiral rebuilding disc 
scenario. 
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Appendix A: Implementation of cooling, feedback 
and star formation in GADGET-2 

To account for cooling, star formation, and feedback processes 
of the ISM, we modified the publicly available version of the 
GADGET-2 code by adding these addi tiona l processes. We f ol- 
lowed the method of lCoxet al] (120061) and ISpringeH d2000h to 
implement these physical process es. We give a brie f description 
below and details are described in lCox et all (120061) . 

Radiative cooling is important for gas to cool down to the 
central region of the dark h alo and then form stars. The cooling 
rate is calculated following iKatz et all dl996l) . in which the pri- 
mordial plasma of H and He and collisional ionisation equilib- 
rium are assumed. Gas particles in these dense regions can form 
stars. The star formation rate is correlated to local gas density 
and anti-correlated to the local dynamic time: 



dp* 
df 



Pgas 



ii 



where 



t d 



'ynamic 



(A.l) 



(A.2) 



The method of deter mining star forma tion rate relies on the 
Kennicutt-Schmidt law (iKennicutti [T998) according to which 

x Pgas' ^ as P^ticles are converted into collisionless par- 
ticles according to the above equation, using a stochastic tech- 
nique. 

To regulate star formation, feedback is needed to prevent gas 
from being consumed by star formation. It assumes that the en- 
ergy released by a supernova is first stored in a new reservoir 
of internal energy. The energy from this new reservoir will pro- 
vide additional pressure to support the gas and prevent it from 
further collapsi ng to form stars. The feedback energy q can also 
be thermalised. lCox et al.l d2006l) introduced two free parameters 
to control this thermalisation process. One gives the time scale 
t feedback of thermalisation process, the other controls the equa- 
tion of state. 

We verified that our implementation of additional routines 
that track the radiative cooling of gas an d star formation. in the 
GADGET-2 code precisely follows that of lCox et al.l d2006l) . We 
have indeed reproduced all the Cox et al figures for all their 
adopted values of feedback (low, medium, and high) and for 
the density dependence of the feedback energy thermalisation 
timescale (n= 0, 1 and 2).. 

We also tested the stability of our mo del galaxies by ru nning 
them in isolation as previously done by ICox et a l. (2006). This 
is shown in Fig. IA.ll which shows how the star and gas surface 
mass density have evolved after 3 Gyr. The gas and star compo- 
nents for both primary and secondary galaxies are stable during 
isolated evolution. This stability is also reflected by the evolution 
of the star formation rate as shown in Fig. IA.2l During the 3 Gyr 
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Fig. A.l. Projected mass density of the primary (left) and sec- 
ondary (right) galaxies for model with mass ratio 3 (model 
M3L34) for which the galaxies are simulated in isolation for ~ 
3 Gyr. The top row shows the stars and the bottom row show the 
gas distribution. A colour bar indicates the surface mass density 
scale shown in the top-left panel. The size of each panel is 160 
by 160 kpc. 
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Fig. A.2. Star formation rates for isolated galaxy models of the 
primary (top curve) and the secondary (bottom curve) for a 
m ass ratio 3 (mode l M3L34). Five times the median feedback 
of ICox etal] (12006) is used in these models. 



evolution, the star formation keeps a roughly constant value. The 
star formation rate of the secondary is about 0. 18 M yr~' and it 
is ~ 1.1 M yr _1 for the primary galaxy. It slowly decreases be- 
cause the star formation affects the gas surface density. The star 
formation rate curve is not smooth because of the hig h feedback 

"etall l2 



value (five times the median feedback of ICox * 



2006) that 



is adopted to regulate the star formation. This test shows that the 
current feedback model can lead to a stable disc and preserve the 
gas before fusion, which is needed to form a disc after a merger. 



Appendix B: Properties of loops formed by 
particles coming back from tidal tails 

lHammer et al.l d2010l) modelled the M31 galaxy with a 3:1, ma- 
jor merger event, which also reproduced the Giant Stream and, 
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for some mode ls, also the northwest bi g loop discovered by the 
PANDA team dRichardson et all 1201 lb . In these models, parti- 
cles returning to the galaxy from the tidal tail can feed the loops 
for several Gyrs ( > 5 Gyr, see their Figure 8). Here we explore 
the properties of such loops, and trace the motion of tidal tail par- 
ticles. In this section, we take an overview of the loop properties 
as they result from our numerical computations. 

Particles stripped into a tidal tail are of two kinds: part of 
them have enough energy to escape the galaxy potential, while 
others reach a maximum distance and then fall back to the 
galaxy. We are mostly concerned by the orbital pattern shape 
of particles coming back from the tidal tail, because they can 
generate many faint features after the fusion. These features are 
very likely observable, because the process c an be maintained 
for several Gyrs (see e.g. lHammer et alll2010l) . These time esti- 
mates are considerably larger than those obtained for encounters 
of two ellipticals and in mu ch better agreement wi th encounters 
of an elliptical and a spiral (Feld mann et al. ll2008l) . This argues 
for the necessity of a cold component in order to produce strong 
and long-lived tidal features. 

A tidal tail formed in a major merger is more complex tha n 
one formed in a minor merger. As shown by IChoi et al . (2007), 
a massive satellite can change the morphology and radial veloc- 
ity of a tidal tail by self-gravity. Meanwhile, the satellite will 
also feel the gravitational attraction from the tidal tail, affect- 
ing its energy and ang ular momentum. The r esonance may af- 
fect the tail properties dD'Onghia et all [20Toh . All these effects 
make tidal tail formation in major mergers more complex than 
in minor mergers and thus more difficult to reproduce. Important 
properties of the orbital pattern of tidal tails include the eccen- 
tricity of loops, which are approximated by ellipses, the width of 
the loops, and the precession rate of the particle orbits. 

In the following we describe the interaction between the 
most massive (the primary) and the less massive (the secondary) 
galaxy in a major merger. To match the NGC 5907 loops better, 
we need to form two tidal tails. In the classic case of an en- 
counter between two prograde discs, we witness the formation 
of two tidal tails, one from each galaxy. On the other hand, in 
our model one galaxy rotates prograde and the other retrograde 
so that both tails form from material stripped from the prograde 
companion, which in our model is the less massive one. The two 
tails form at different times, namely the times of the first and 
the second closest passages. The motion of all the particles that 
constitute the tidal tail is set before the interaction by the orbit 
of the secondary and its disc spin. After being stripped from the 
secondary by the gravitational interaction, these particles keep 
some orbital information from their initial orbits, i.e. from their 
progenitors. When the secondary is close to pericentre where the 
tidal perturbation is maximum, the tidal force will add e nergy to 
the particles of the secondary ( Springel & Whit el 1 1 9991) , whic h 
makes its particles unbound and ejected (D'Onghia et al., 2010). 
Alon g the tidal tail, par t icle binding energy increases with radius 
(see iHibbard &MihosL fl995t iDubinski et all [l996), i.e., parti- 
cles ejected further from the centre have higher energy, leading 
to an increasing fallback time. 

The motion of particles falling back from the tidal tail can 
be easily understood since the dynamical friction force can be 
neglected due to their low mass. If we also neglect the time evo- 
lution of the central potential, the motion is simply that of test 
particles in a static central potential, which is a rosette (Binney 
& Tremaine 1987, hereafter BT87). The angle between two con- 
secutive apocentra depends on the potential, being degrees in a 
Kepler potential and 180 degrees in a harmonic potential. Thus 
the orbit can be approximated by a precessing ellipse, while the 



precession angle can give information on the potential of the 
merger remnant. 

Figure lBTl shows an example of the particle locations at dif- 
ferent epochs. The tidal tail particles have been tagged at 3.5 Gyr 
after the beginning of the simulation, i.e. after the fusion of the 
two galaxies (see Sec. [5}, so that we can follow their individual 
motions. These particles continue to feed the loops as they return 
from the tidal tail. At 4.2 Gyr the first loop is already formed, and 
then at 5.6 Gyr a second loop is added. As time increases further, 
the third and fourth loops are formed, as shown by the 8.6 Gyr 
panel of Fig. lBTTl 

Particles are continuously feeding the loops following a sim- 
ilar track from the tidal tail. They first enter into the first loop and 
then continue by orbiting in higher order loops. Because parti- 
cles falling back later have higher energies (or elevations rela- 
tively to the mass centre), the size of the first loop increases with 
time, as well as that of other loops after their formation (see Fig. 
IB. II ). We also fitted the first and second loops with ellipses at 5.6 
Gyr and 8.6 Gyr as shown in Fig. IB. 21 Owing to the permanent 
motion of particles, it is not always simple to separate particles 
from different loops. We used the distance of particles from the 
centre at the tail-forming time (3.5 Gyr panel of Fig. IB. U to iso- 
late particles for each loop. In Fig. |R2] different colours indicate 
different distances to the centre. The green particles are initially, 
at the time of forming the tail, closer to the centre and now form 
the second loop, while the black particles that are farther from 
the centre form the first loop. This shows that the particles form- 
ing the different loops come from different distances to the centre 
at the time the tidal tail is formed. In this simulation we notice 
that the first two loops formed by particles returning after reach- 
ing their apocentra are approximated well by two ellipses with a 
phase angle ~ 80 degree difference due to precession, as shown 
in the left hand panel of Fig. IB. 31 

The loop width depends on the orbital properties of the tidal 
tail. It has been shown that resonance is important for form- 
ing a long tidal tail for a spin system. Strong resonance oc- 
curs when rotation and orbital frequencies obey Q^, tr = Q. or bu 
dToomre & Toomrel 1 1972k [Barnes & HernquisA fl992l) . As long 
as we can control the tidal tail length and width, we may control 
the shapes of the loops at some level. For various merger geome- 
tries, we can decrease the loop eccentricities by decreasing the 
orbital eccentricity of the secondary, because they are related. 
This is illustrated in the right hand panel of Fig. IB. 31 where we 
compare the loop eccentricity with the orbital eccentricity of the 
second-order loop. Loops used in this figure are coming from 
models with different mass ratio, pericentre, eccentricity, bary- 
onic fraction (from 6% to 9%), and gas fraction. The arrow in- 
dicates the region for which orbital eccentricities provide loop 
eccentricities consistent with the observations. 



14 



Jianling Wang et al.: Loops formed by tidal tails as fossil records of a major merger 



3.5Gyr 


4.2Gyr 




S;6G^- 


9.5Gy^- ... 








tT! ■ . -V^ 1 ■ 





Fig. B.l. Loops formed by the particles falling back from second tidal tail. Each panel is centred on the centre of mass of the 
secondary galaxy, indicated by a cross. In this figure, the tidal tail is formed from a 3:1 merger (model M3L34, see section[5]). The 
size of each panel is 300 by 30 kpc. The stellar mass surface density of the loops at 8.6 Gyr is consistent with the estimation of 
iMartm ez-Delga do et alJ (120081) for NGC 5907. The top row shows the loop formation process as seen from the observed viewing 
angle, while the bottom row gives a face-on view of the loops at the same times. The red arrows in the second column indicate 
the first loop, the pink arrows in the third column show the second loop, the blue and green arrows in the fourth column show the 
third and fourth loops, and the black arrow in the last column shows the fifth loop. In the bottom-left panel, solid and dotted circles 
indicate the radii of the original location of the particles which constitute the tidal tails at 5.6 Gyr and 8.6 Gyr, respectively. Particles 
with small elevations are selected from the region between the green and the black circles, and those with high elevation, from 
beyond the black circle. 
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Fig. B.2. First and second loops for the same model as in Fig. IB. II are fitted by ellipses for simulation with mass ratio 3. The left 
panel is at 5.6 Gyr and the right panel shows loops at 8.6 Gyr. Particles belonging to different loops are naturally selected by their 
distance to the centre in the tail-forming epoch (3.5 Gyr). Green indicates particles with small elevations (see Fig. IB. II bottom-left 
panel), and black points indicate large elevation particles. Green points populate the second loop, while black points show particles 
of the first loop. Particles coming back from the tidal tail begin to outline the first loop and then enter in the second loop. The 
parameters of the fitted ellipses are given in the top left corner of each panel. They are the semi-major axis a, semi-minor axis b, 
focus distance, eccentricity, and position angle. 
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Jianling Wang et al.: Loops formed by tidal tails as fossil records of a major merger 
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Fig. B.3. The left panel gives the distribution of loop phase angle differences between first and second loops. The right panel 
compares the loop and orbit eccentricities. Solid black line is the one-to-one relation. This statistics is coming from various models 
with different mass ratio (3, 4, and 5), orbital geometry and dark matter fraction. 
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